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Abstract 



We discuss vector, axial-vector, scalar and pseudoscalar two-point func- 
tions at low and intermediate energies. We first review what is known from 
chiral perturbation theory, as well as from a heat kernel expansion within 
the context of the extended Nambu-Jona-Lasinio (ENJL) model of ref. [p!^ ]. 
In this work we derive then these two-point functions to all orders in the 
momenta and to leading order in l/Nc- We find an improved high-energy 
behaviour and a general way of parametrizing them that shows relations be- 
tween some of the two-point functions, which are also valid in the presence of 
gluonic interactions. The similarity between the shape of the experimentally 
known spectral functions and the ones we derive, is greatly improved with re- 
spect to those predicted by the usual constituent quark like models. We also 
obtain the scalar mass Ms = 2Mq independent of the regularization scheme. 
In the end, we calculate fully an example of a nonleptonic matrix element in 
the ENJL-model, the vr"'' — ir^ electromagnetic mass difference and find good 
agreement with the measured value. 



1 INTRODUCTION AND REVIEW OF KNOWN RESULTS IN QCD 



We shall be concerned with two-point functions of the vector, axial-vector, scalar 
and pseudoscalar quark currents: 

V^\x) ^ q{x)^,^q{x) (1) 

4")(x) = q{x)^^^,—q{x) (2) 

Si-){x) ^ -q{x)^q{x) (3) 

P^^\x) ^ q{x)i^,^q[x) (4) 

where A^"^ are Gell-Mann SU(3)-matrices acting on the flavour triplets of Dirac 
spinors: q = {u{x), d{x), s{x)). Summation over the colour degrees of freedom of the 
quark fields is understood. These are the quark currents of the QCD Lagrangian 
with three light flavours u, d, s, in the presence of external vector Vf^{x), axial vector 
afj,{x), scalar s{x) and pseudoscalar p{x) field matrix sources; i.e.. 



with 



and 



Cqcd{x) = Cqcd + + 750^)? - - ^Tsp)? , (5) 



^QCD = E G^^G^^^^^ + iq^id, + iG,)q (6) 

^ A=l 



G,^9sT. ^Gi^\^) (7) 

A=l ^ 

the gluon field matrix in the fundamental SU{Nc — 'i)coiour representation, with 
G^^^ the gluon field strength tensor 

G'-^^ix) = - - gsfABcG'f ^G^f^ , (8) 

and Qs the colour coupling constant (ct^ = gl/Ai:). 
We want to consider the set of two-point functions: 



(9) 
(10) 

(11) 



KA<l)ab = i J 


1 d^xe"^-"" 


<0|t(\/W(,x)\/W(0)) 


|0 > 


^iu{Q)ab = i J 


[ d^xe^^-^ 


< o\T {a(;\x)a^^\o)) 


|0 > 


K(^)ab = ^ J 




< 0\T (v^''\x)S^''\0)) 


|0 > 
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nj'(g)„b = i J < 0|T (A(f)(a;)p(^)(0)) |0 > (12) 

n^iq)ab = I J c/'xe'"-" < 0|T [S^''\x)S^^\0)) |0 > (13) 
U^{q)ab = I J d^xe'"-^ < 0\T (^P^''\x)P^^\0)) \0 > . (14) 

The relevance of these two-point functions to hadronic physics, within the frame- 
work of current algebra, was emphasized a long time ago (See refs.[|l| to 0].). With 
the advent of QCD it became possible to compute the short-distance behavior of 
these two point functions using perturbation theory in the colour coupling con- 
stant because of the property of asymptotic freedom: as(Q^ = — >> Aq^j-,) ~ 
log~^{Q'^ / Aqqjj), (See refs.|^] to [^.). Non-perturbative corrections at large ap- 
pear as inverse powers in 0. The inclusion of these corrections, combined with 
a phenomenological ansatz for the corresponding hadronic spectral functions at low 
energies, has developped into the active field of QCD sum rules. (See e.g. ref.@ for 
a review and references therein.) The QCD behaviour of the two point functions 
above at very low values, is controlled by chiral perturbation theory (xPT). It 
is important for our later purpose to review here what is known in this respect. 

From Lorentz-covariance and SU(3) invariance, the two-point functions above 
can be decomposed in invariant functions as follows: 

KMat = iq,qu - q'g,.M\Q')Sab + q.quI^P iQ')Sab (15) 
Ki^U = iq.Qu - q'g,uM\Q')5ab + q,qun^A\Q')Sab (16) 



{(l)ab 
{(l)ab 

Il'{q)ab 

n^(g)a6 



Ili{Q')q,Sab 
n^(Q2)ig/,(5„6 

Ils{Q^)Sab 

HpiQ^ab . 



(17) 
(18) 
(19) 
(20) 



The low energy behaviour of these invariant functions in SU{2)l x SU(2)^ xPT has 



been worked out in ref. |T^. It is easy to extend their analysis to SU{3)l x SU{3)ii 
XPT. In the chiral limit, where the quark mass matrix Ai 0; and with neglect 



of chiral loops, which are non leading in the l/Ai'c-expansion [|T^, the results are as 
follows: 



n?^(Q^) 



(21) 



(22) 
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= ^ - 4(2iJi - Lio) + 0{Q') (23) 
n?^(g^) = (24) 
n^(g^) = (25) 

n^(Q^) = ^ (26) 

n^(Q2) = 85^(2^8 + H2) + O(Q') (27) 

n^(g2) = ?||/| + sBl{-2L, + i/s) + 0(g') . (28) 
Q 

The functions H^-*, n£f and 11^ get their leading behaviour from the pseudoscalar 
Goldstone pole. The residue of the pole is proportional to the pion decay constant 
/o (/o ~ /,r = 93.3MeV^). The constant Bq is related to the vacuum expectation 
value 

< 0|gg|0 >\,=uAs= -foBo (1 + OiM)) . (29) 
The constants Lg, Lio, Hi and H2 are coupling constants of the O(p^) effective 



chiral Lagrangian in the notation of Gasser and Leutwyler |12[. The constants 
Lg and Liq are known from the comparison between xPT and low energy hadron 
phenomenology. At the scale of the p mass: 

Lg = (0.9±0.3) X 10"^ (30) 

and 

Lio = (-5.5±0.7) X 10"^ . (31) 

The constant Hi and H2 correspond to couplings which involve external source fields 
only and therefore cannot be extracted from experiment unambiguously. 

It has been recently shown that the extended Nambu Jona-Lasinio model 
(ENJL-model) describes the values of the low energy parameters rather well. In 
ref. |T^, various relations between the parameters of a low energy effective La- 
grangian were obtained that were independant of the input parameters and possible 
low energy gluonic corrections. Among these was the relation 
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f^M^ - flMl = fl (32) 

between couplings and masses of the lowest vector and axial-vector mesons and the 
pion-decay constant. This corresponds to what is usually called the first Weinberg 
sum rule |l|]. The second Weinberg sum rule, in the lowest resonance saturation 
form i.e., fyMy = f\M\, was, however, not satisfied exactly though the deviation 
was numerically small. Part of the motivation which triggered our interest in the 
present work has been to understand the origin of this unsatisfactory result. As 



we shall show, the low energy expansion method used in ref.|13] is inappropriate to 
draw conclusions about the intermediate Q^-behaviour of two-point functions. The 
relevant contributions can however be easily summed using a Feynman-diagram 
expansion of the four fermion couplings in the ENJL-model without introducing 
collective fields. Section 3 gives the details of this summation method. The same 
results could be obtained of course using the collective field method, provided though 
that all higher orders in the Q^-expansion were kept. 

It is well known that the Weinberg sum rules play a crucial role in the calculation 
of the electromagnetic 7r+-7r° mass difference in the chiral limit. In a previous cal- 



culation of this mass difference [14| within the framework of the effective action ap- 
proach [jl5l , it was shown that the quality of the matching between long-distance and 
short-distance contributions to the photon-loop integral is still rather poor when 
the vector and axial-vector spectral functions are replaced by constituent quark 
spectral functions alone. We shall (see also ref. [0) discuss this problem again in 
section 4, within the framework of the ENJL-model and with the full Q"^ /Mq depen- 
dence summed. This calculation, which as we shall see is rather succesful, provides 
the first non-trivial example of a genuine one loop calculation in the ENJL-model. 
There are other applications one can now think of doing; in particular calculations 
at the next to leading order in the l/A'c'-expansion. We plan to investigate this in 
the future. 

Throughout this paper we shall work in the chiral limit. The inclusion of cor- 
rections, whenever neccesary, due to nonzero current quark masses can however be 
done with the present technology. 

The rest of the paper is organized as follows. In section 2 we give an overview of 
the ENJL-model and the two-point functions in this model in the low-energy limit. 
We also describe here the parametrization usually used to go beyond the low-energy 
expansion in this model. Section 3 is the mainstay of this work. Here we derive 
the two-point functions to all orders in the momenta. We also discuss how gluonic 
corrections can be included and present numerical results. In the next section we use 
these two-point functions to calculate completely in the ENJL-model a non-leptonic 



4 



matrix-element, the 7r+ — vr*^ electromagnetic mass difference. In the last section we 
present our main conclusions. In the appendix we derive the underlying relations 
that allow for the same type of results for the two-point functions, as those found 
in ref. for the low-energy constants. 



2 TWO-POINT FUNCTIONS IN THE ENJL-MODEL 



2.1 A brief review of the ENJL— model 



The scenario suggested in ref. |jT3[, assumes that at intermediate energies below or 
of the order of the spontaneous chiral symmetry breaking scale A^^,, the Lagrangian 
of the ENJL-model is a good effective realization of the standard QCD Lagrangian 
i.e., ^ 

^^QCD ^QCV + ^^N^L + ^'n'jL + ^(tI 

X 



(33) 



with 



and 



5,P 
N.JL 



[x] 



iVcA2 



<^ X a.b 



(34) 



(35) 



Here a and b are u, d, s flavour indices and colour summation within each quark 
bilinear bracket is implicit; = |(1— 75)g(x) and g^j, = |(H-75)g(x). The couplings 
Gs and Gy are dimensionless quantities. In principle they are calculable functions 
of the ratio of the cut-off scale A^ to the renormalization scale Aj^g. In practice the 
calculation requires knowledge of the non-perturbative behaviour of QCD and Gs 
and Gy will be taken as independent unknown constants. In choosing the forms (0) 
and (|3^) of these four quark operators we have only kept those couplings that are 
allowed by the symmetries of the original QCD Lagrangian, and which are leading in 
the l/A^^c-expansion |Tl|. With one inverse power of Nc pulled-out, both couplings 
Gs and Gy are 0{1) in the large Nc limit. The A^ index in C^^d means that only 
the low frequency modes of the quark and gluon fields are to be included. 

The basic assumption in considering the ENJL-model as a good effective La- 
grangian of QCD is that at intermediate energies below or of the order of the spon- 
taneous chiral symmetry breaking scale, the operators and >C]vjl are the lead- 
ing operators of higher dimension which, due to the growing of their couplings Gs 
and Gy as the ultraviolet cut-off approaches its critical value from above, become 
relevant or marginal. 

As is well known in the Nambu Jona-Lasinio model jl^, the C^jj^ operator, for 
values of Gs > 1, is at the origin of the the spontaneous chiral symmetry breaking. 
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It is this operator which generates a constituent chiral quark mass term [U is a 
unitary 3x3 matrix which collects the pseudoscalar Goldstone field modes): 



- MgiqLU^qn + qRUqL) = -MqQQ , (36) 
like the one which appears in the Georgi-Manohar model JTSf; as well as in the 



effective action approach of ref. ||15|| 



As discussed in ref. [0, the C^j^ operator is at the origin of an effective axial 
coupling of the constituent quark fields Q{x) with the Goldstone modes 



with 



IsAQr-i'^^Q , (37) 



and r(0,e) the incomplete gamma function: 

r(0,e)=/ -e-^. (39) 

For e small, 

r(0,e) = -loge-7i? + O(e) , (40) 
with = 0.5772 ■ ■■, Euler's constant. We recall that in eq. 



= I {i^ [d, - i{v, + a,)] i-i[d,- i{v, - a,)] e) ; (41) 

with 

U = ii (42) 

and 

Q = Ql + Qr (43) 

Ql = iqi , Qr = ^qR ■ (44) 

The appearance of incomplete gamma functions is due to the proper time regu- 
larization which is used in calculating the fermion determinant via the Seeley-De 
Witt expansion. The type of integrals which appear are (e = Mg/A^) 



^ ^;r^r(n-2,e) 



167r2 (M3)"-2 
n= 1,2,3,- •• (45) 
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The result for the axial coupling in eq. pq ) is the result to leading order in 
the l/A^^c-expansion. In terms of Feynman diagrams it can be understood as the 
infinite sum of constituent quark bubbles shown in Fig. la, where the cross at the end 
represents the pion field. These are the diagrams generated by the Gy - four fermion 



coupling to leading order in the l/A'^c-expansion [|T9]. The quark propagators in 
Fig. la are constituent quark propagators, solutions of the Schwinger-Dyson equation 
in the large Nq approximation, which diagrammatically is represented in Fig. lb. 

In its simplest version where one assumes furthermore that all the relevant gluonic 
effects for low energy physics can be absorbed in the couplings Gs and Gy, the 
ENJL-model has only three free parameters: Gs, Gy and A^^. We find it useful 
to specify these in terms of Mq, and qa instead. For this purpose one should 
remember that (e = Mq/ 



er(-l,e) =e-^-er(0,e) 



(46) 



2.2 Low Q^— behaviour of two— point functions in the EN JL— model. 

As we discussed in the previous section, the low-energy behaviour of two point 
function quark currents, is governed by the constants Bq, Lg, Liq, Hi and H2. 



In the ENJL-model, these constants have been calculated in ref. with the results 

p2 ^'O' . ,^2 



AMi^gATiO,e) 




2Hi-L 



167r2 8 



(47) 
(48) 
(49) 
(50) 

(51) 



In fact, the calculations made in ref.[T^ suggest a possible improvement of the 

) when the contribution from the corresponding 



invariant fuctions in eqs. (^) to ( 
resonance propagators is also taken into account, with the results 



n?^(g^ 



-A{2Hi + L 
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nS^(g') = T^-4(2iJi-Lio)- 



M2 + Q2 ' 



(52) 
(53) 
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and 



ml 



2Lq + H2 



(54) 



For the other invariant functions, the results are the same as in eqs. ( pTD to (Pq) 
with the parameter values given in eqs. ( ^7|) to (^Ij). Several comments on these 
results are in order: 

"""""" n . rR-rn 



i) Two relations which follow from the ENJL-model [O 



2H, - Lio = -fl/2 and 2Hi + L,o = -/y/2 



(55) 



These relations are of the same type as the first Weinberg sum rule relation in 
eq. (p2D i.e., they are independent of the input parameters and possible low energy 



gluonic corrections. As first shown in ref . , they are crucial to ensure that the low 
energy effective theory is compatible with the known short-distance properties of 
the underlying theory - QCD. It is reassuring that the ENJL-model indeed respects 



these constraints. Using these relations, we can also write Uy^ and 11^'' in the form 



n 



n 



(1) 

A 



(1) 
V 



+ 



(56) 



(57) 



a form much more similar to the usual vector meson dominance (VMD) phenomeno- 
logical parametrizations found in the literature. 

ii) The constant in eq. (p3|) denotes the coupling 



(58) 



in the effective scalar Lagrangian. Here 



X = 2Bo [s{x) + ip{x)] 



(59) 



with s{x) and p{x) the external scalar and pseudoscalar matrix field sources. As 
discussed in ref.|jl3|, the couplings Lg and H2 receive contributions both from the 
quark loop - which we denote Lg and H2 - and from the integration of scalar fields 
- which we denote Lf and iff; i.e.. 



-^8 + Lo 



and 



Ho = Ho + HS . 



(60) 



In fact 
2d 



Ml 



2LI + Ho = 2Ho 



Nc 1 



_ . r(o,e) 

167r2 4^^r(-l,e)2 ^ ^ 



■he) 



21(0, e)]^ 



(61) 
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iii) Equations (|52|), ( p3D and ( p4D imply a specific form of tfie resummation to all 
orders in an expansion in powers of . As we shall see, this is not however the 
correct form which follows from the exact resummation of Feynman diagrams, to 
leading order in the l/A'c'-expansion. 

We shall finally give the expressions for the masses M^, My and Ma which are 



obtained in the ENJL-model in ref.[13 



3 r(0,e) 

Ml = QMl— . (64) 

^ r(o,e) 

3 FULL g2_j3EpENj3ENr[, TWO-POINT FUNCTIONS IN THE 
ENJL-MODEL 

3.1 The Vector Two— Point Function 

To illustrate the method, we shall first discuss with quite a lot of detail the vector 
invariant function T{y\Q'^). In the ENJL-model, and to leading order in the I/Nq- 
expansion, we have to sum over the infinite class of bubble diagrams shown in Fig. 
2a. Algebraically, this corresponds to the sum 

{q^q'' - q'gn^l^y^ + (g^?" - q^gH^^ (^]^) ^ 2(g«g^ - q^g^W^v^ + ■■■, (65) 

where the explicit factor of 2 in the second term comes from the two possible con- 
tractions between the fermion fields of the vector four-quark operator. The overall 
result IS 

Notice that for this two-point function, only the vector four quark interaction with 
coupling Gv can contribute. No mixing between different operators can occur in this 
case. The one loop bubble in Fig. 2b corresponds to the bare fermion-loop diagram 
of the mean field approximation defined in ref. |T^, which in what follows we shall 
denote with an overlined expression 

'U^^\q') = U^^\q') . (67) 

It is easy to see that at the n-loop bubble level, the corresponding expression for 
^ni,^^(Q2) will be given by the n — 1 - loop bubble result multiplied by the coupling 
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Gv and one more factor of the one-loop result, i.e., 



'li'v'm = -'K\Qr-^i'Q')K'iQ') ■ (68) 

X 



This series can then be summed with the result 



V - ^ m • (69) 

g^^n«(g2) 



We discuss next the calculation of IIy\q^) in some detail. The spectral function 
associated to Ily\Q'^) is the one corresponding to the QQ intermediate state in a 
P-wave and can be calculated unambiguously with the well known result 



i/™II<;>W^^i l.^Wl-^«(<-4MS). (70) 



167r2 3 V t J ^ t 

The function Ily\Q'^) we seek for has to obey three criteria: 

i) It must obey the relevant Ward identities. 

ii) Its discontinuity should coincide with the spectral function in eq. (|70D. 

iii) When expanded in powers of it must reproduce the heat kernel calculation 
of the effective action approach, with the same proper time regularization results. 
These three criteria will in fact apply to all the two-point functions we discuss. 

To proceed with the calculation oiIly\Q'^), we write a once subtracted dispersion 
relation for this function 

K'm = n^\o) - r T7T?^-^-n«(t) . (71) 

Jo t t + Q"^ TC 

Requirement iii) implies that Ily^O) is fixed, with the result (e = Mg/A^) 

n?\o) = -f^|r(o,6) . (72) 

With ^Imlly\t) in eq. ([70|) inserted in the integrand of the r.h.s. of eq. (|7lD ; and 
with the successive change of variables 

— ^ = 1 - and ?/ = 1 - 2x , (73) 

we have that 

dt g2 I 



16.2 - 2-)^ [1 + 2x(l - X)] ^.^^t(l-x) • 
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In order to match with the proper time regularization which has been used in the 
calculation of ny''(0), we next replace the denominator in the r.h.s. of eq. ( [f4D as 
follows 

I , ^^p-r[A/2+Q2a.(l-x)] /^c-N 



Performing an integration by parts in the x-variable we then finally get the result 



/o t t + Q^n' 

where xq is a short-hand notation, which we shall use from here onwards, for 

_ + Q^x{l - x) 



(77) 



Combining eqs. (|7lD, (|72|) and ([76|), we obtain 

Tlv iQ') = ^8£dxx{l-x)T{0,XQ) . (78) 
The first few terms in a Q^-expansion of this expression are 

in agreement with the proper time regularized heat kernel effective action result. 

The imaginary part of Ily\Q'^) evaluated from eq. (|78D using the ie prescription 
Xq (m^ + Q'^x{l - x) - iej /A^ in the log term: T{0,xq) = -logxg - 7s + 
0{xq), reproduces the spectral function in eq. (|70D. 

We shall now try to cast the result for Ily\Q'^) in eq. (^) in the simple VMD- 
form of eq. (^) : 



where we have set 



and 



V 



^ AGy Jq dxx{l — x)T{0, Xq) ^ ^ 

We find that the full Q^-dependent vector two-point function can indeed be cast 
in the VMD-form of eq. (|56D provided that the coupling parameters fviQ"^) and 
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MviQ"^) become dependent. Their value at = happen to coincide, in this 
case, with the couphngs in the low energy effective Lagrangian i.e., 

MJW=^0)^MJ^?|^. (83) 



and 



m' = ^) = fv = ^\me). (84) 



The product fy{Q'^)My{Q'^) is scale-invariant. 

In order to see the hadronic content of the full vector two-point function we 
propose to examine the complete spectral function 



TT 



2 



t^Jmn«(t) 



n 2 



^5) 



+ 

and plot it as a function of t for the input parameter values 

Mq = 265MeV, = llQhMeV (86) 

and 

QA = 0.61 . (87) 



These are the values corresponding to fit 7^1 in ref. ||13|]. The plot is the one shown 
in Fig. 3 (the full line). For the sake of comparison, we have also ploted in the same 
figure the spectral function ^Imliy\t) corresponding to the mean field approxima- 
tion (the dashed line). The improvement towards a reasonable simulation of the 
well known experimental shape of the = 1^ , 1 = 1 hadronic spectral function is 
rather notorious. 

3.2 The Axial- Vector Two-Point Function 

The infinite series of bubble diagrams we have to sum in this case is formally 
very similar to the one already discussed in the previous subsection. Again, only 
the vector four-quark interaction with coupling Gy contributes in this case with the 
result 

n\"(Q^ 



1 + Q2 80.n^i)(g2) 



The axial two-point function 11^ (g) in the mean field approximation, has however 
more structure than the corresponding vector two-point function because now, due 
to the presence of a constituent quark mass, the axial invariant function 11^^ in the 
decomposition corresponding to eq. (|T6]) doesn't vanish. Both 11^ (Q^) and 11^ (Q^) 
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have associated spectral functions which can be calculated unambiguously with the 
results' 



l,nm)-^'-f^i-'-fe,-ur,) (90) 

Notice that the three spectral functions ^Imlly\t), ^Imll^^\t) and ^Imll^^\t) 
satisfy the identity: 

-Imlty\t) - -/mn^'\t) - -ImTifit) = . (91) 
71 n n 

In fact this is nothing but a particular case of a general Ward identity which two- 
point functions in the mean field approximation must obey: 

Tl^\Q')-Tl^l\Q')-uf{Q') = 0. (92) 

The proof of this identity can be found in the Appendix. It is precisely this identity 
which guarantees that the first Weinberg sum rule in the mean field approximation 
is automatically satisfied. 

From the asympotic behaviour of Imll^^\t) in eq. (pO| ) we conclude that the 
dispersive part of the function n^'*(Q^) obeys an unsubtracted dispersion relation. 
To this we have to add the pole term calculated in the effective action approach i.e., 

tt(o)/^2n -2/^ , /-"^ dt 1 (0) 

^AiQ) = ^ + l ^-^-n^W' (93) 



with (e = Mil Is?' 



/^ = ^4M^r(0,e) . (94) 



Using the same change of variables as in eqs. (|73|) , we obtain 

Next, we use the same proper time representation for the denominator in the right 
hand side as in eq. (|75D , and perform an integration by parts in the x-variable, with 
the result 

with xq defined in eq. (|77D. Inserting this result in the r.h.s. of eq. (p3D leads to the 
final result 
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which defines a running /^(Q^) in the mean field approximation. The first few terms 
in a Q^-expansion of /^(Q^) are 



fliQ') = ^4M5 |r(0, e) - lr(l, e)^ + Om \ , (98) 



in agreement with the proper time regularized heat kernel effective action result. 
Once we have calculated II^^Q^) and Ily^Q^), the function II^^Q^) follows 



from the Ward identity in eq. (|9^) with the result 



^a\Q') = I dxT{Q,XQ) + 1^ dxx{l - x)r(0,XQ) [ . (99) 



■(1), 



It is nevertheless instructive to calculate 11^ {Q ) independently as an illustration of 
the method we are using. First we observe that the dispersive part of n^''(Q^) needs 
a subtraction. On the other hand, we know from the effective action calculation that 
Tf^{Q'^) has a pole term (see eqs. (Bl) and (pi)) 



We recognize the second term in the r.h.s. of this expression as the constant term in 
the Q^-expansion of /^(Q^) in eq. (P5|), which means that the subtraction constant 
needed for the dispersion relation is only the third term. However this term is 
precisely the same as the one corresponding to the vector function ny''(0) in eq. (|7^). 
We must therefore separate the spectral function ^ImU.^^ [t] in two pieces: one 
which reproduces the vector spectral function ^Iml\.y\t), for which we shall write 
a once subtracted dispersion and the rest. But this is precisely the Ward identity 
separation we already pointed out in eq. (^). We then have 




^ r:^i^(i + ziii£iwi_l^fi0(t_4M2). (101) 

the two dispersive integrals are those calculated before in eq. (p6|) for the unsub- 



tracted piece. Putting these results together leads to the result in eq. (p9|). 

With these results in hand, we shall now try to cast 11^'' (Q^) in eq. (|88|) as close 
as possible to the VMD-form of eq. (^T]): 



tt(1) tt(0) 
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From the calculated expression for ufiQ^) in eq. (P7D, it follows that 

^ - Q'^0^^a\Q') = l + 4Gy^ I'dxmxQ) . (103) 

At = 0, the r.h.s. is precisely g^^ (see eq. (|38|)). This is not a surprising result. 
The evaluation of the axial vector form factor of a constituent chiral quark from the 
infinite series of bubble graphs in Fig. la, leads to the result 



1 

l + {Gv/AlHMi^Io'dxT{0,XQ) 



which at = coincides with the axial coupling constant qa obtained in the 
calculation of the low energy effective action in ref. |TB[. Using this result, we can 



now rewrite the r.h.s. of eq. ( |102| ) in the following simple form: 



where 



MTmtV ^ (105) 



m') = 9A{Q')fl{Q') ; (106) 

MliQ') = — (Q2) . (107) 



and 



with /^(Q^), M^{Q^) and f^{Q^) as given in eqs. {^), (||) and (|8l|), respectively. 
Notice that at = 



= 0) = f: = ^4Mj^^^r(0, e) ; (109) 



i.e., the same expression which appears in the low energy effective action; however 
MKQ^ = 0) and f\{Q^ = 0) do not coincide with the expressions for the couplings 
Ml and f\ as given in eqs. (|6^) and (^51), (^9]). We now have instead 

M1(Q2 = 0) = — M2 and fl{Q' = 0) = g\f^ . (110) 

gA 

The remarkable new result is that now, in terms of the running couplings and running 
masses, both the first and second Weinberg sum rule are satisfied 

fUQ')MUQ') = fl{Q')Ml{Q') + fl{Q') (111) 

and 

fUQ')Mt.{Q') = fl{Q^)M\{Q^) . (112) 



15 



The fact that f\{Q'^ = 0) does not coincide with the f\ defined in the effec- 
tive action approach can be understood from a comparison between eqs. (|57|) and 
eqs. (|105|) . Both expressions coincide to O^Q"^): 

n?le,. (i) = ^ + 2/l + W) (113) 

n«|e,. ® = ^ - i^^^^r(i, e) + 2/l(g2 = o) + o{q') (ii4) 

It is the fact that part of the constant term is reabsorbed in the Q ^-dependence of 
that is at the origin of this difference. The precise relation is 



/l = /l(Q' = 0)-^^<7ir(l,e) (115) 



3.3 The Scalar Two-Point Function 

The full scalar invariant function from the sum of the infinite series of bubble dia- 
grams has the form 

It only involves the four-quark operator with Gg-coupling. As in the previous vector 
and axial-vector discussion, we now proceed to the calculation of the scalar two- 
point function Ils{Q^) in the mean field approximation. The associated spectral 
function can be calculated unambiguously, with the result 



From its assymptotic behaviour at large t, it follows that the function Ils{Q'^) obeys a 
dispersion relation with two-subtractions for the term proportional to tJl — AMq/t 



and one subtraction for the term Mq^I — AMq/t. Accordingly, we write the dis- 
persion relation 

T[s{Q') = n5(o) + (n'^(o) - ^^r(i, e] 



:ii8) 



t2 t + Q2 IQ^2 

where 115(0) and 11^(0) are already known from the effective action calculation i.e., 

n5(0) = ^4MJ [r(-l, e) - 2r(0, e)] ; (119) 
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and 

n;(o) = -^2(r(o,6)-^r(i,e)) . (120) 

Notice that only the divergent pieces (i.e., terms proportional to r(— 1, e) and r(0, e); 
but not r{n,e), ~ 1) are retained in the subtraction constant. The integral we 
have to compute is 



dt 1 / 4M|, 1 /"i , ,o 1 , , 

TTTWr - ^ = 2 i ''^(^ - mITqmT^ ■ ''''' 

where we have made the standard change of variables of eq. (ffSl) . Using the rep- 
resentation of eq. (|75D, and doing one integration by parts in the variable x, and 
replacing eqs. (|119|) and ( |120| ), we obtain the result (e = xq = — — 



X X 



n.(g^) = ^4M2|r(-i,e)- 2(^ + 1) f^dx^x^)^ . (122) 

We shall next evaluate the full IlsiQ'^) function in eq. ( p,16| ), and try to cast it in 
a form as close as possible to eq. (|5^): 



NcK 



2-] 



Using the fact that Gg^ = er(— l,e), we can cancel the terms proportional to 
r(— l,e) in the denominator; and write Ils{Q^) in the simple form 

UsiQ') = SBl I2Z3 + k + ' (124) 



where 



and [xq = .2 ) 

X 



M's = 4M^ , (125) 



' Ml - iQ.-^^^2s^,dxmxQ)- 2Minm- ^ ^ 

The term 

%Bl{2U + H,) = -^4M5r(-l, e) = ^ (127) 

is not quite the same as in the effective action calculation. The value of Ils{Q^ = 0) 
however, coincides with the same result as in the effective action: 

Us{0) = ^Bl{2L, + H,) = ^4MJ^^;^ [r(-l, e) - 2r(0, e)] . (128) 
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As in the previous subsection, it was an advantage to reabsorb some of the constant 
terms into the pole term, compared with the expression obtained from the effective 
action calculation, to obtain a simple expression to all orders in Q^. 

The striking new feature of the summed scalar propagator is that the scalar mass 
is constant and Ms = 2Mq. This is to be contrasted with results of previous work 
in the literature, see e.g. ref. and references therein. Eqs. (|125|) and (|127|) are 
also true in the presence of gluonic corrections. They, and the fact that Ms = 2Mq, 
are a consequence of the identities derived in the appendix. It is only the specific 
form of the functions that depends on the inclusion of gluons or not. 



3.4 Two— Point Functions with Mixing 

The case of the other two-point functions n^'*(Q^), Il^{Q'^) and Il^{Q'^) is 
somewhat more involved because they mix through the Nambu-Jona-Lasinio four- 
fermion interaction terms in eqs. (|3^ ) and (pS]). Therefore, for this case the result 
at the n bubble level is a matrix equation in terms of the two-point functions at the 
n-1 bubble level: 



fgvTif 9sK 
















X 




/ 


""'Hp J 



(129) 



where Ily\Q'^), n^^((5^) and Ilp{Q^) denote, as usual the two-point functions at the 



one loop level calculated in the mean field approximation; and gv, gs are short-hand 
notation for 



9v 



NcAl 



Q and gs 



NcAl 



The series we have to sum, in a vector-like notation, is 

oo oo 



n=l 



n=l 



(130) 



(131) 



and the 



where B[Q'^) denotes the 3x3 two-point function matrix in eq. 
n((5^)'s three component two-point function vectors. This series can be summed. 



l-B{Q^j 

The matrix {1 — B) can be inverted, and after some algebra we get the result 



A(Q2) 



[1 - gsU'')uf + gsiKf 



A(g2 



(132) 

(133) 
(134) 
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and 



np(Q') 



{l-gvIif)Iip + gv{jfM) 



(135) 



with A((5^) the function 

A(Q2) = (^1 _ ^^nj^) (l - (75np) - gs9viTil?j' . (136) 

It is illustrative to see what the result would be for the scalar two-point function 
n'^(Q^) if we had carried the analysis keeping the functions Yiy^ and H^: 

{l-gyUf){l-gsTls)+9sgv{Ki) 

If we now set Ily'* = = 0, we recover the result of eq. ( |116|) . 

We need now to calculate np(Q^) and lij^iQ'^). To calculate lip^Q"^) we proceed 
as for the other two-point functions we have already calculated. The corresponding 
spectral function is 



l^rnYlAt) = §,SMi'^ 1 - '-^e(t - m-) . (138) 



For the large t-behaviour we conclude that Ilp{Q'^) obeys a dispersion relation 
subtracted twice: 

Up{Q') = np(0) + Q'U'AO) + Q' / T^T-7^-ImUp{t) , (139) 

with the divergent pieces (divergent when e = Mg/A^ — > 0) of llp{0) and np(0) as 
known from the effective action calculation i.e., 

np(0) = ^4M2r(-l,6) ; (140) 

and 

nU0) = --f%2r(0,e). (141) 

The integral in the r.h.s. of eq. (|139|) has already been calculated in the section 
about the scalar two-point function. We then find the result {xn = ^q+Q ^(-^ 

X 

np(Q') = ^{4Mjr(-l,e)-2g2^'rfa;r(0,XQ)} . (142) 

Let us check that this result satisfies the three criteria we discussed in section 
3.1. First, there is a Ward identity which relates Ilp{Q'^) to 11^^ (Q^), because of 
the axial current divergence condition 

d^^iQYl^Q) = 2MqQi^,Q (143) 
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in the mean effective field theory. The Ward identity in question, which we proof in 
the appendix, is 



AM^Up{Q')+AMq < QQ >= (Q^yufiQ^) ; 



(144) 



(145) 



and [|T^ 

<QQ >=-^4Mjr(-i,6). 

Equations (p7| ) and ( |142| ) indeed satisfy this Ward identity. Second, the spectral 
function calculated from eq. ( |142|) via the ie prescription: 



- tx(l -x)-ie 



log I 



+ me 



' txil -x)-Ml 



is the same as in eq. ( |138| ). Finally, the first few terms in the Q^-expansion of 
np(Q^) coincide with those calculated in the effective action approach with the 
proper time heat kernel regularization. 

The last two-point function to calculate is I\.^.j{Q'^). The corresponding spectral 
function is 



1 _p A^^ 

TT loTT^ 



4M^ 



-Oit - 4M2 ) . 



(146) 



The function Il]^j(Q'^) obeys a once subtracted dispersion relation. We already 
have encountered the same situation with one of the terms in the scalar two-point 
function, which we have discussed in detail. Therefore, we give the final result only 



N, 



c 



167r2 



4Mq / dxT{0,XQ) . 



(147) 



We can now proceed to the explicit calculation of the functions 11^^ (Q^), 11^^ (Q^) 
and np(g2) in eqs. (|13|), ([13|) and (|13|). We find that 



l-gsHpiQ')] uf{Q')+gs{uUQ')y = 



and hence 



n?)(Q^ 







(148) 



(149) 



a result which must hold in QCD at the chiral limit. Equation ( |148| ) also implies 
that 

Q' 1 



AiQ') = 1 - gsUpiQ' 



and therefore 



n^(g^) 



M2 2r(-l,e) 



< QQ > 



dxT{0, xq) 



Q' 



(150) 
(151) 
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As for Ilp{Q^), from eq. (|135| ) and using the relations (|148[ ) and (|150| ) above, we 
find 

"""2)- i-9,np«j^) -M^ + yf(QW- 

We see from this resuh that mixing between the pseudoscalar and the longitudi- 
nal axial degrees of freedom occurs at all orders in the Q^-expansion. Because of 
eq. (|150| ), Ilp{Q^) has now a pole at = 0. We want to stress the fact that the 



final full results in eqs. (|149|) , ( |151j ) and ( |152| ) are very different to those at the 



one-loop level approximation. It is only when all the contribution to leading order 
in 1/Nc are summed that these relations, which are expected features of QCD, ap- 
pear. These results are a big improvement with respect to the QCD effective action 
approach at the mean field approximation. 

3.5 Inclusion of Gluonic Corrections 



In Ref. |]T3[ the low-energy corrections due to the lowest dimensional gluonic con- 
densate were also explicitly included. A general analysis based on the possible types 
of terms here corresponds essentially to keeping all the overlined functions as un- 
determined parameters but satisfying the relations derived in the appendix. This 
follows from the fact that gluonic lines connecting different fermion loops in Fig. 2a 
are suppressed by extra factors of l/Nc compared to the leading contribution. 

The correction due to the leading gluonic vacuum expectation values can in fact 
be easily included by using the results for two-point functions calculated for use in 
QCD sum rules P]. These corrections can be rewritten in terms of the dimensionless 
parameter 

9 = 7^i-G') (153) 



and the set of functions 



"1 1 
Jn= I dx — . (154) 



The corrections needed for the spin-1 parts are (for Nc = 3) 

27r2Q4 



27r2Q4 

For the scalar two-point function we need 



ns(Q') = (-1 - 2J-1 + 3 J2) . (156) 



47r2Q4 
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The spin-0 axial-vector, pseudo-scalar and mixed two-point functions have very large 
cancellations in the denominator and are numerically very unstable when gluonic 
corrections are included. They can be handled similarly in principle. The gluonic 
correction terms do also satisfy the relations derived in the appendix as required. 

We have checked that using the above formulas the two-point functions including 
non-zero gluonic vacuum expectation values converge for small values of to the 
low energy expansion with these corrections included; i.e. those of eqs. (21)- (28) 
with the values of the parameters calculated including gluonic corrections. All the 
nice features of the two-point functions as given in the previous subsections are still 
valid since the underlying cause for these properties were the relations derived in 
the appendix and the gluonic corrections have to satisfy those as well. 



3.6 Numerical Results 

In this section we plot the two-point functions as calculated in the previous sub- 



sections. As described in subsection p.5| we can also include gluonic corrections. In 
view of the result of ref. [1^ that a very good fit to the low-energy parameters could 
also be obtained without gluonic corrections, we only show the effect of the gluonic 
corrections in the vector two-point functions. The input values used for all of the 



plots are those of fit 1 in ref. 0. They are Mq = 0.265 GeV, = 1.165 GeV 
and = 0.61. The variation with the input parameters can be judged from table 
2 in ref. [l^. The size of the changes here is similar to the ones obtained there at 



values of momentum transfer = 0. 

In fig. 4 we have plotted the vector-two-point function for positive values of 
for \^Cp from to 1.5 GeV. The full line corresponds to 11^ (g^) in the full ENJL- 
model. The dashed line is the corresponding result using the effective approximation, 
eq. (|56[). The vector meson mass for the values of the parameters Mq, A^ and qa 
which we have fixed, is about 0.81 GeV. The short-dashed line corresponds to the 
vector-two-point function in the QCD effective action model of ref. [|T^]. As can 
be seen in the figure, the full resummation leads to lower values for the two point 
functions than those of the low-energy formulas when extended to higher Q^. In 
the same figure 4 we also show the effect of the gluonic corrections. The dotted line 
uses the same input values as given above but has a non-zero value for the gluonic 
background, we have set g = 0.5 (see eq. ( |153| )). It can be seen that the effect of 
gluonic corrections is large at small values of but grows smaller at higher values 
of 

In fig. 5 we have plotted similarly to the vector case, the result for Q^U^{Q^). 
(The extra factor of is included to remove the pole at = due to pion ex- 
change.) In contrast to the vector two-point function, there is also a significant 
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difference between the one-loop result of eq. ( pTj) (dashed line) and the full resum- 
mation of eq. (105) at low Q^(the full line). This is due to what is usually called 
the pseudoscalar-axial- vector mixing and was in our previous work described by the 
coupling qa- The value of this two-point function at = determines The 
dashed hne corresponds to the two-point function using the effective approximation 
of eq. (|57|) . The axial- vector mass here is about 1.3 GeV . 

Fig. 6 shows how the summation of the whole series of diagrams has produced the 
pole at = that is required by the spontaneous breaking of the chiral symmetry 
for the pseudoscalar two-point function. The one-loop result (dashed line) does not 
have this behaviour, but the full resummed version (full line) of eq. ( |152| ) does. The 
short-dashed line is the low-energy extrapolation of eq. (28). Here we see how the 
full resummation correctly reproduces the low-energy behaviour as derived in ref. 
13, but for larger values of momentum transfer, it starts to differ appreciably. 



We have not plotted the scalar two-point function. The pole is, as we have proven 
above, always at Ms = 2Mq. This pole is generated by the full resummation. 
Neither have we plotted the mixed pseudoscalar-axial-vector two-point function 



since this has the very simple behaviour of Eq. ( 151 ). 

The overall picture of the high energy behaviour of the ENJL-model that emerges 
after the resummation, is improved compared to the behaviour obtained from the 
low-energy expansion. This is illustrated by the fact that now it also satisfies the 
second Weinberg sum rule. The ENJL-model has another advantage over the simple 
QCD effective action model of ref. |T5|. By virtue of the extra 4-quark interactions 
present, this model naturally contains more or less correct meson poles while the 
simple quark version, that corresponds to the one-loop result (or essentially the 
use of the overlined two-point functions), does not. This means that for positive 
values of the two-point functions are considerably enhanced both in the real and 
imaginary parts as compared to the one loop result. The importance of this type of 
behaviour can be seen, e.g. in the determination of some low-energy constants using 
dispersion relations. As an example we show in fig. 3 how the imaginary part of the 
vector two-point function gets enhanced considerably over the one-loop result. 

The advantage of using the full ENJL-model over a parametrization with meson 
resonances is that the number of free parameters remains within limits. For instance, 
if one tries to extend the analysis of ref. [l^ to non-leptonic matrix elements using 
a parametrization with vector mesons it requires the knowledge of weak decays of 
vector mesons which have not been observed experimentally. 

In ref. |l^ the masses and couplings of the mesons were determined from the low- 
energy expansion. These are essentially given by various combinations of derivatives 
of the two-point functions at = 0. An alternative way of determining the meson 
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masses is to determine them by looking for the poles and residues of the full two- 
point functions. This procedure can be questioned on the grounds that for euclidean 
momenta quark confinement is not so important but for momenta of > 4:Mq we 
get effects of free quarks included. The two-point functions still have poles though. 
As an example we give the position of the poles for the vector, axial-vector and 
scalar for the parameters used above, in table |l]. The pion mass is of course exactly 
zero in both cases since we work in the chiral limit. In general the masses are lower 
than those derived from the low energy approximation. 

4 THE 7r+-7r° ELECTROMAGNETIC MASS DIFFERENCE 

To the lowest order in the chiral expansion, the effect of virtual electromagnetic 
interactions to lowest order in the fine structure constant = e^/47r, generates a 



term in the effective action without derivatives 22 



d'^x [e^CitrQU{x)QU^{x)} , (157) 

where U is the unitary matrix which collects the pseudoscalar Goldstone fields and 
Q the quark electric charge matrix Q = l/3diag(2, —1, —1). Expanding eq. ( |157|) in 
powers of pseudoscalar fields 



2e^i 

J n 



e'CitrQUQU^ = -^^(tt+tt^ + K+R-) + 0{4>^) , (158 



one sees explicitly that this term leads to a 7r+-7r° (and K^-K^) mass splitting 

J TT 

The constant Ci, like is not fixed by symmetry requirement alone. It is deter- 
mined by the dynamics of the underlying theory. Formally, it is given by the integral 
representation 



^e'C, = -ze^ / (g^g- - gV^)nS(g^) , (160) 

where 

na = i(ni;'-n<") , (lei) 

with Uy^ and H^"* the vector and axial vector invariant functions we have discussed. 
Performing a Wick rotation in eq. (|160| ) leads to the sum rule 
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The purpose of this section is to discuss the evaluation of the Am^ sum rule 
above based on the two-point function results discussed in the previous sections. It 
is then convenient to split the Q^-integral into long-distance {0 < < /i^) and 
short-distance (/x^ < < oo) parts: 



dQ'--= dQ'---+ Q'-- (163) 

Jo JtJ.^ 



We shall concentrate first on the long-distance part calculation. 

a) Long— distance contribution. Phenomenological approach 

The very low contribution to the integral 

is fixed by chiral perturbation theory (see eqs. (0) and (p3D): 

u\^\q') - n«(Q2) = ^ - 8L,o + 0{Q') , (165) 



from which it follows that [14 




fi^ + 0{fx^)} . (166) 



The known correction term 0{fi'^) in the parenthesis of the r.h.s. can be used to 
estimate the value of the /z^-scale at which we can trust the validity of the xPT- 



contribution. From the fact that 



^ ~ (167) 

we conclude that the xPT-result in eq. ( |166|) can only represent correctly the long- 
distance contribution to Am^ up to scales 

iApt < Ml . (168) 

Obviously, this is too small a scale to trust numerically a direct matching with the 
short-distance contribution, which as we shall see later, it is expected to be valid for 
/i^-scales larger than a few GeV^ at least. (See however the first paper in ref. P2 
andref. ig.) 

Since the early work of Das et. al. 0, the traditional phenomenological approach 
to the calculation of (Am^)^^^ has been to include the effect of vector and axial- 
vector particle states in the Q^-integral, using a parametrization that is constrained 
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to satisfy the first and second Weinberg sum rules. The usual phenomenological 
VMD-model parametrization is 

- (169) 

and 



n;^^ = ^ + , (170) 



^(1) _ 2^ 2/iMi 

g2 ^ Ml + Q2 

with the constants /y, /I, My and M\ constrained by the relations 



fl + flMl = flMl (171) 

and 

PyMt. = flM^, (172) 
which ensure the convergence of the limits 

lim Q2(n« _ ny ) ^ and lim Q\U^y^ - H^^^) ^ ; (173) 

i.e., the super convergence relations which lead to the first and second Weinberg sum 
rules. One then has 

For Ma, My ^ oo, with /i^ fixed we recover the first term of the xPT calculation 
in eq. ( |166D . If we let the scale /x^ go to infinity; then, for M^ = \/2My, one finds 
the early result of Das et al. : 

(Am^)g = ^lMpog2 = 1.4 x lO^MeV^^ . (175) 

Experimentally, 

(m^+ - m^o)^^p = (4.5936 ± 0.0005)Me1/ , (176) 
while the phenomenological result of Das et al. corresponds to 

(m^+ - m^o)g = 5.2MeK , (177) 
Recent phenomenological evaluations of the Am^ sum rule, which include explicit 



chiral symmetry breaking effects, can be found in refs.|23] to [25 
b) Long— distance contribution in the ENJL— model 



The calculation of (Am^)^^ in the QCD effective action approach of ref. [p!5[ 



which corresponds to the mean field approximation of the Nambu Jona-Lasinio 
model, was reported in ref. Jl^. It is the approximation where 

n« - n"' ^ n« _ n»' = ,,r(o, , (i78) 

26 



which is the result obtained in eqs. ( p2D and (^). This leads to the result (e 



X 



r(o,6 



In ref. |Tj| a proper time regularization for the photon propagator was used; and 
for simplicity, the /i^-scale was identified with A^. The shape of this mean field 
approximation evaluation versus /i^, for the input value of Mq and which we have 
been considering (i.e., the value corresponding to fit 1 in ref.|]l3|: Mq = 2Q5MeV 
and A^ = 1165Me\/) is plotted in Fig. 7. 

The evaluation of {AmDiD in the full ENJL-model, with the expressions of the 
two-point functions IIy\Q'^) and II^^\Q'^) obtained in the previous section leads to 
the result 

with the Q2_dependent functions M^{Q^), M\{Q'^) and fl{Q^) as given by eqs. (jS^), 
(ITOTD , ([104|) and ([106|) . 

The shape of {/S.m'^)ENjL versus /i^ is shown in Fig. 7. We expect the integrand 
in eq. ( |18U| ) to be a good representation of the low and intermediate energy scales; 
and therefore, the matching with short-distance evaluation should now be much 
smoother than in the case of the mean field approximation. This we discuss in the 
next subsection. 

c) Short— distance contribution and numerical results. 

In QCD perturbation theory I{y\Q'^) = 11^'' (Q^) . Spontaneous symmetry 
breaking induces a deviation from this result which, at large and to leading order 
in the l/Ai'c-expansion, can be calculated using the operator product expansion, 



with the result (M and first ref. in 22 



(^<;'-^»') = -iJfM!)««.»^ (181, 

where {Nc — > 00): 



ge 2 TT 
NcasiQ') 



:i82) 



and 

< ^7/>(Q2) >= < >(log(Q/AQCD))* . (183) 

Inserting this asymptotic estimate in the short-distance expression for Am^, leads 
to the result 
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Fig.7 also shows the shape of {Am'^)sD versus fi^ for various values of the invariant 
quark condensate < i/jiIj >. Obviously, as the scale becomes small (Am^)^^ 
diverges. The matching between (Am^)^^ and (Am^)^^ is defined by the optimal 
choice of fi^ which minimizes the variation of the total Am^. As seen in fig. 8 this 
occurs at value fi ^ 950 MeV; and in fact around the value, the stability is rather 
good. The corresponding value of Am^ in this range, is 

Ami ~ 1-3 • 10"^ GeV^ , (185) 
and agree well with the experimental value, the horizontal dashed line in fig. 8. 

5 CONCLUSIONS 

In this paper we have extended the general analysis of the ENJL-model as done 
in ref. [jl3| beyond the low-energy expansion. We have calculated directly the two- 
point functions within the ENJL-model to all orders in momenta. The relations 
that the one-loop results have to satisfy lead after the full resummation to a set of 
rather simple forms for the two-point functions. It should be stressed once more 
that these are satisfied independent of the gluonic interactions and are thus valid in 
a wide class of ENJL-like models. 

The resulting expressions are, for the vector-axial- vector cases, very similar to the 
ones usually obtained assuming some kind of vector, axial- vector meson dominance. 
The full resummations have a well behaved high-energy behaviour. They satisfy 
both the first and the second Weinberg sum rules. The resummation also obeys the 
Ward identities of the full theory. 

Simple expressions were also found for the other two-point functions. A byproduct 
was a proof that within this class of models the scalar two-point function always has 
a pole corresponding to a mass of twice the constituent quark mass. Our derivation 
only depends on the underlying symmetry properties of the Lagrangian and is hence 
regularization scheme independent. The full resummation also reproduced the pole 
at = in the pseudo-scalar two-point function explicitly showing how this model 
obeys the Goldstone theorem. 

Finally, the two-point functions derived were used to start evaluating nonleptonic 
matrix elements within the class of ENJL-like models. We have estimated the elec- 
tromagnetic 7r+ — 7r° mass difference and found good agreement with the measured 
value. 

*Thc continuous curve is the one corresponding to the choice | < tpi/j > \ — (281Mey)'^, which is 
the value predicted in the ENJL-model for the input values Mq — 265MeV and = 1162Mey. 
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APPENDIX 

In this appendix we derive the Ward identities that the one-loop two-point func- 
tions have to satisfy. We first give a derivation based on the heat-kernel expansion 
and a general analysis of the type of terms that can contribute to the two-point func- 
tions. This method allows for explicit contact to be made with the regularization 
chosen in the heat-kernel expansion. A second method is essentially the traditional 
way of deriving Ward identities but we have to take into account that {qq) ^ 0. 
The second method can also be used to derive some of the identities that the full 
two-point functions have to satisfy. 

The one-loop two-point functions are calculated using the Lagrangian ([/ = 1) 

L — qilj> q — Mqqq — q{s — ipj5)q — qVq . (186) 

The last equality is the definition of V and the covariant derivative contains 
the vector and axial-vector external fields. The real part of the effective action in 
Euclidean space using the heat kernel expansion is then given by (e = Mq/A^): 

^'f^ ^ " 3^^ 5/^'' " ^' / d'^'^^nix) . (187) 

The Hn{x) are the Seeley-DeWitt coefficients and these are constructed out of 
E, R^^ and their covariant derivatives. These are defined by 

V^Ve = -D^D^ + E + M^ and [D^, Dj\ = R^, . (188) 

In terms of the external fields s, p, and they are (only terms that can contribute 
to two-point functions are given): 

E = i-fi^'jsMQ (r^ - Z^) - -a^^uRnu 

+ Mqs + p^ + -f^d^s - , (189) 

Rtxv = --^{l^iv + r^^-'^z{l^^-rf,^)) . (190) 

Here we see that E and R^^, vanish for vanishing external fields so only terms con- 
taining at most two factors of E and R^i, can contribute to the two-point functions. 
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The first two coefScients are: 

Hq = 1 and ni = -E . (191) 

Til thus contributes to the scalar and pseudoscalar two-point function. These are 
the only two-point functions that contain a quadratic divergence. 

The 7i„>2 only contain two types of terms that can contribute to two-point func- 
tions. Let us look at all possibilities. 

Terms with a single E. These are of the form D'^^"'~^^E and are total derivatives, 
so they do not contribute to the two-point functions. The same argument applies 
to terms with a single Rf^^. 

Terms with one E and one P^jy. Extra derivatives acting on these can always 
be commuted, the commutator introduces an extra factor of R^u and then only 
contributes earliest to a three point function. We can also use partial integration. 
All this type of terms can thus be brought into the form 

D,,D,ED'^^-^^R^, = i [D^, D,] ED^^^-'^R^, . (192) 

The commutator becomes an extra factor of R^i, so this type of terms does not 
contribute to two-point functions. We conclude that the TCn>2 only contribute to 
two-point functions through terms like 

^^2{n-2)^ ^ DaRai3D^^'"~^^Df,R^p and P„^D2("~2)i?„/3 • (193) 

Using Eq. ( |19CI| ) the last two terms are of the form 

D^VapD^^'^-'^D^v^p + D^a^f.D^^'^-'^D^a^p (194) 

and 

VafsD^^^'^^v^^ + a^fsD^^'^-^^a^f, . (195) 

So these contribute only to the transverse part and equally for the vector and the 
axial-vector two-point function. The first term has a part containing R^^^y as well. It 
contributes only to the transverse part, and equally for the vector and axial- vector 
case. 

The remaining type of terms can be rewritten using the explicit form of E. 

+16MjP92("-i)p + im^S + . (196) 

The axial-vector terms contribute only proportionally to g^j^p. This together with 
the above contribution leads to: 

Tt^^\Q') = nf^)(Q2) = 0, (197) 

n?\Q^) = nT{Q')+nf{Q'). (198) 
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The first two of tliese equations appear because tlie vector current in tlie Lagrangian 
Eq. ( |186| ) is conserved. Tfie tliird one is tlie reason wliy tlie first Weinberg sum rule 
is satisfied even at the one-loop level. It also guarantees both Weinberg sum rules 
after the resummation. Including the contributions from Tii we also have 

-2MqUI,{Q') = QmfiQ'), (199) 
2Mq'Up{Q^) = -2(QQ)-Q2nf,(Q2) ^ (200) 

UsiQ') = Up{Q')+Qmf{Q') . (201) 

In Eq. (|200|) we have used the relation between the coefficient of Tii and the quark 
vacuum expectation value. In the chiral limit this vacuum expectation value is 
determined uniquely by the contribution of Tii = —E. This derivation is also valid 
in the presence of low-frequency gluons. The effective action after including the 
low-energy gluonic effects through gluonic vacuum expectation values, still has to 
be constructed out of E and R^^. This was precisely the argument used in Ref. [jl3| 



to obtain relations between the low-energy coupling constants that are independent 
of the gluonic corrections. The results following from the relations ( |197[ - p01| ) after 
resummation are the equivalent relations for the two-point functions. This is what 
we used in Sect. 3 to rewrite all the two-point functions in terms of essentially two 
functions and one constant. 

The preceding derivation was obtained using the Seeley-DeWitt expansion to all 
orders. Let us now show how several results can also be obtained from the underlying 
relations in the Lagrangian (|186|). These relations are: 



d^Qlf^Q) = 0, (202) 

df, (g7M75g) = 2iMQqj5q , (203) 

{CH^),?S(0)} = d'^X^d^x). (204) 

Eq. ( p04| ) is vahd at equal times. a,b are colour-flavour indices and a,j3 are Dirac 
spinor indices. 
We start from 

g.nj: = Jd'x{d,e^''-){0\T{A,{x)A^{0))\0) (205) 



-2Mq J d*x{0\T{P^{x)A,{0))\0) 

- I d^a;5^o5(x°)(0| [A^(x), A,(0)] |0) (206) 



= -2iMQl]!f'^ . (207) 

The matrix element of the equal time commutator vanishes for two identical currents. 
This follows from Eq. ( |204| ) . Putting in the form of the two-point functions this 
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leads to 

q\uf = 2MQqJlP^^ (208) 



or the same as equation ( p.99| ). In the full theory we have d/^^A^ = so the identical 
derivation leads to: 

uf{Q^) = 0. (209) 

This equation is satisfied by the fully resummed two-point function. 
A similar derivation leads to 

g^nj = 2zMq!F- J dV5^o5(x°)(0| [A^{x),P{0)] |0) . (210) 

Here the equal time commutator worked out using Eq. ( p04| ) does not vanish. A 
term proportional to the quark vacuum expectation value remains and leads to Eq. 
( [2001 ). In the full theory df^A^ = so we obtain 



nfM)(Q^) = -2^. (211) 

This equation is also satisfied by the fully resummed two-point function. 
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Table 1: Values of the masses determined from the poles in the two-point functions 
and from the low-energy expansion of ref. . 



Meson 


ref. 


13 


1 


Pole 


My 

Ma 
Ms 


0.81 GeV 
1.3 GeV 
0.62 GeV 


0.70 GeV 
0.9 GeV^ 
0.53 GeV = 2Mq 



' In the resummed version there is an strong enhancement around this value of the 
two-point function. It does not become a pole with the values of the parameters 
chosen here. 



Figure Captions 

Fig. 1: (a): The set of diagrams summed to obtain gA^Q^)- X is the insertion 

of the pion field and the other lines are fermions. 
(b): The gap equation. The thick line is the full fermion propagator. 

The thin line is the bare fermion propagator. 
Fig. 2: (a): The set of diagrams to be summed for the two-point functions. 

(b): The one loop fermion bubble. 
Fig. 3: The spectral function of the vector two-point function. The full line is the 

full ENJL result. The dashed line is the result at one- loop. Here Q = 
Fig. 4: The real part of the vector two-point function. Plotted are the full result 

with (labelled gluon) and without (labelled full) gluonic corrections. The 

VMD parametrization (eff) and the one loop result (1-loop). 
Fig. 5: The real part of the axial- vector two-point function multiplied by Q^. The 

effect of axial-pseudoscalar mixing that the full resummation reproduces is 

visible at all Q's. The labels have the same meaning as in Fig. 4. 
Fig. 6: The real part of the pseudoscalar two-point function multiplied by Q^. No- 
tice how the resummed version produces the pole at Q = 0. The labels have 

the same meaning as in Fig. 4. 
Fig. 7: Curves for Am^ in terms of the scale fi. Plotted are the long-distance 

part for the pion exchange term only (LD-CHPT), the result of ref. 

1^ (LD-mean) and the ENJL- result after the resummation (LD-ENJL). 



The short distance contributions are plotted for three values of {QQ) = 

-(194 Me1/)3(SD194), -(220 MeVf and -(281 MeVf. 
Fig. 8: The full result Am^ versus fi corresponding to the sum of the long distance 

ENJL result with the short distance evaluation using (QQ) as given by the 

ENJL-model (see text). 
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